Quantum fluctuations of a vortex in a dilute Bose-Einstein condensate 
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A vortex in a quasi two-dimensional Bose-Einstein con- 
densate is subject to the Magnus force and can be effectively 
described as a planar particle in a uniform magnetic field. 
Quantization of this effective particle leads to the lowest Lan- 
dau level where the most localized wave function is a gaussian. 
In this gaussian state vortex position seems to fluctuate with 
an average magnitude set by the magnetic width of the gaus- 
sian. We readdress this problem using the number-conserving 
version of the Bogoliubov theory. We find that the Bogoliubov 
mode that might be interpreted as a fluctuation of vortex po- 
sition actually does not contribute to position fluctuation at 
all. The only non-zero contribution comes from phonons but 
it is an order of magnitude less than the simple estimate, 
based on the magnetic width of the effective gaussian wave 
packet. 

PACS: 03.75.Fi, 05.30.Jp, 32.80.Pj 



I. INTRODUCTION 

Bose-Einstein condensation of dilute alkali gases [1] 
opens up new possibilities to investigate Bose-Einstein- 
condensed systems under flexible and well controlled con- 
ditions. For example it is possible to manipulate with the 
dimensionality of the trap to achieve quasi one or two- 
dimensional condensates [2]. Temperatures can be re- 
duced to the pK range [3] bringing the system close to the 
quantum regime. Various manifestations of superfluidity 
have been subject of intensive research, the most spec- 
tacular have been probably vortices [4] and Abrikosov 
vortex lattices [5] in rotating traps. It is predicted [6] 
that further increase in the rotation speed of a quasi-two- 
dimensional trap can drive the system from the regime of 
BEC to the regime of strongly correlated bosonic quan- 
tum Hall liquids with even denominator filling fractions. 

A dilute bosonic gas in a rotating two-dimensional trap 
at very low temperature is a setup required for quantum 
Hall liquids but at the same time it can be used to ex- 
perimentally test quantum fluctuations of a 2D vortex in 
a dilute Bose-Einstein condensate. In fact melting of the 
Abrikosov vortex lattice driven by quantum fluctuations 
of vortices is the first step from a Bose-Einstein conden- 
sate towards strongly correlated quantum Hall liquids. 
In this paper we develop a theory of quantum vortex 
fluctuations. 



II. VORTEX FLUCTUATIONS 

An effective action for a planar vortex has been worked 
out in Ref. [7]. In terms of vortex position on the plane 
X = (X,y) the action reads 



S cS = J dt ^X + ^hp e z (X x 1) 



(1) 



Here m e s is an effective mass of the vortex that depends 
on the size of the system, p is planar density of atoms in a 
uniform condensate, the integer k is the winding number 
of the vortex, and e z is a unit vector along the z-axis 
normal to the plane. Formally this is the action of a 
planar particle in a uniform magnetic field normal to the 
plane, as we can see from its Euler-Lagrange equation: 

rricfiX = nTip X x e z . The effective Lorenz force on the 
right hand side is the Magnus force. Because vortices 
with |k| > 1 are unstable, in the following we consider 
only k = +1. 

The effective particle can be quantized and its eigen- 
states can be classified in the degenerate Landau lev- 
els. The energy gap between Landau levels is so 
in the limit of vanishing effective mass the spectrum 
can be truncated to the lowest Landau level (LLL). The 
degenerate eigenmodes of the LLL can be chosen as 
(Z*Y er^ pZ * z . Here Z = X + i y is a, complex po- 
sition of the vortex. We get the most localized eigenstate 
when the non-negative integer 2 = 0: 



e -yz-z = e 



i P (x 2 +y 2 ) 



(2) 



This gaussian seems to determine the quantum fluctua- 
tion of the vortex position. 

In the next Section we rederive this result in the frame- 
work of the particle-number-conserving Bogoliubov the- 
ory [8], but then we use the same Bogoliubov theory to 
show that this result is incomplete and even misleading. 



III. VORTEX FLUCTUATIONS 
IN THE BOGOLIUBOV THEORY 

In the mean-field approximation a stationary conden- 
sate wave function is a solution of the time-independent 
Gross-Pitaevskii equation 



^^-—V^ + Ngo^U + Vcj) . (3) 
Am 
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Here go is an effective 2D interaction strength and N 
is a total, conserved number of atoms in the box. In 
this section we ignore the trap potential V and consider 
a uniform 2D condensate in a box. In the strong in- 
teraction Thomas-Fermi (TF) regime (Ngo = J > 1) 
the uniform condensate has a wave function <j) = \J p/N 
with p, = pgo ■ In this case vortex has a wave func- 
tion (j>o(r, 9) = f(r) exp(i6) in polar coordinates r and 0, 
where the modulus f(r) interpolates between /(0) = 
and /(oo) = yfp~/N. 

Small fluctuations around the vortex can be expanded 
into Bogoliubov modes with complex amplitudes b m 



= N 

where 

Af-Hb, b*) 



h(x) + J2 b mUm(x) + b* m V* m {x) 



(4) 



1 + J2m,n b* m bn(u m \u n ) + b m b n {v* n \u n ) + 



b* m b* n {u m \vl) + b m b* n {v n \v m ) (5) 



is the normalization factor, such that (</>b|</>b) = 1. In ad- 
dition to the usual phonon modes with the characteristic 
phonon dispersion relation, there is a normalizable mode 
with zero frequency 




uo = -\ — 



V n = - 
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(6) 



Here z = x + iy = re tS . This mode is normalized in the 
usual way (uo\uq) — (vo\vo) = 1. When truncated only to 
the zero mode the Bogoliubov expansion (4) becomes 
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(7) 



This expression can be interpreted as the leading terms 
of a Taylor expansion in powers of the complex vortex 
position Z = X + iy = ^/|" b . To sec quantum fluc- 
tuations of Z we must go to the quantum version of the 
iV-conserving Bogoliubov theory. 

Any A^-particle state can be written as a quantum 
superposition over condensates with different condensate 
wave functions 4>b 



-I 



d 2 b^p(b,b*) \N:<j) b } 



(8) 



Here d 2 b is an abbreviation for FJ m d 2 b m . The state \N : 
4>b) is a Fock state with N atoms in the wave function (f>b 



from Eq.(4). For N » 1 the action of the quasiparticle 
annihilation operator 



a v. 



N 



(9) 



with u m = {u m \i})) and v m ■ 
turns out to be very simple: 

b m \ip) = 



(Vmli'), on the state \ip), 



d 2 b^{b,b*) b m \N-.fo) 
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\N : h) ■ (10) 



The annihilation (creation) operator bm acts on a 
wave function ip(b,b*) as a differential operator. A l b- 
representation' of the operators is 



i 1 (v. d 

^ = V^r^dbfr 

f t _L( h * d 

m V2 V m dbrr 



(11) 

(12) 



It is easy to check, that [&m,frjj = 5 mn . 

The Bogoliubov vacuum state |0(,) is the state with- 
out any Bogoliubov quasiparticles, so its wave function 
is annihilated by all the operators b m : 



b m \%) = 0, 



it implies: 



-6*6 



(13) 



(14) 



If b m is applied l m times to the vacuum state in or- 
der to create l m quasiparticles, then one gets the l m - 
quasiparticle wave function (^/2b m ) Lm e~ b b . It is quite 
remarkable that these wave functions look just like the 
orbitals of the lowest Landau level. This may be a lit- 
tle bit surprising because after all, they are the eigen- 
states of the quasiparticle harmonic oscillators and not 
of any particles in a magnetic field. Also, more generally, 
a quasiparticle Fock state with definite numbers of dif- 
ferent quasiparticles has a characteristic 'polynomial x 
exponential' form (JJ m {b* m ) lm ) e~ b b . A general wave 
function is tp(b, b*) — P(b*)e~ b b with an arbitrary poly- 
nomial in 6*'s. 

As we already could see in Eq.(7), in the case of our 
planar vortex one of the fictitious particle coordinates 

'bo- 



bo materializes as a complex vortex position Z = 

The Bogoliubov vacuum wave function (14) for this mode 
becomes 



(15) 
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and it is identical with the wave function (2). In contrast 
to Eq.(2) the wave function (15) is derived from first 
principles without any reference to the effective action 
(!)• 

In the following we use the same first principles to show 
that Eqs.(2) and (15) are misleading. It docs not make 
sense to truncate the theory to the vortex zero mode 
because this mode is not independent from other Bogoli- 
ubov modes. 

IV. A CONTRADICTION 

The quantum fluctuations described by Eq.(15) result 
in quantum depletion of atoms from the condensate wave 
function 4>o(x). The density of atoms depleted from the 
condensate in the state (8) with the wave function (15) 
truncated to the vortex mode can be worked out as 

dpo(x) I \u (x)\ 2 + \v (x)\ 2 . (16) 

Here uq and vq come in a symmetric way like in Eq.(4). 
This dpo(x) 7^ at r = 0, the non-zero contribution is 
coming from the term |uo(x)| 2 . The fluctuating vortex 
position results in a non-zero density of depleted atoms 
inside vortex core. This is quite natural: a vortex with a 
fixed position has an empty core but when the position 
is uncertain then the average density of atoms inside the 
core in non-zero. There is a one to one correspondence 
between the magnitude of position fluctuations and the 
density of depleted atoms in the core. 

On second thought, Eq.(16) appears strange. After 
all the density of depleted atoms in \0b) state, using the 
Bogoliubov theory, is given by the sum X)m=o l w m(^)| 2 
which contains all v m } s but no w m 's. When truncated 
only to the vortex mode the sum becomes 

dp (x) = M£)| 2 (17) 

in apparent contradiction with Eq.(16). Both Eq.(16) 
and Eq.(17) result from a truncation of the Bogoliubov 
theory to the vortex mode. However, the two truncations 
are done in a different way. The resulting contradiction 
shows that such a truncation does not make sense. 

The reason why it makes no sense is that different Bo- 
goliubov modes (u m ,v m ) are not independent from each 
other. Their wave functions are not mutually orthogonal 
in the usual sense: for example, in general (u m \u n ) =/= 
for m ^ n. The wave function of the vacuum (14) factor- 
izes implying independence of different 6 m 's, but the non- 
trivial correlations between different Bogoliubov modes 
are hidden in the TV-particle condensate \N : ^) making 
the mode independence rather elusive. 



V. PHONONS DRIVE VORTEX 
FLUCTUATIONS 

In contrast to Eq.(16) equation (17) shows vanish- 
ing density of depleted atoms at r = 0, i.e. there are 
no vortex position fluctuations originating from the vor- 
tex mode. At the same time the contradiction between 
Eqs.(16) and (17) shows that we cannot judge on posi- 
tion fluctuations using any of the two truncations, but 
we have to use the full expression 

oo 

dp(r = 0) = ]T K(r = 0)| 2 (18) 

m=0 

including all Bogoliubov modes. We know that the vortex 
mode (to = 0) has no contribution to dp(r = 0). To 
work out the contribution from the other modes in a 
controlled way we put the system in a 2D harmonic trap 
with a frequency u>. In the dimcnsionless trap units the 
stationary Gross-Pitaevskii equation becomes 

M 0=-iv 2 0+ir 2 + . 9 H 2 0. (19) 

Its vortex solutions </>o = J{r)e 10 for several different val- 
ues of g are shown in trap units in Fig. 1 . 
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FIG. 1. Profiles of the vortex wave function for several 
values of the interaction strength g. The solutions were ob- 
tained by relaxation to the ground state with vortex symme- 
try. Except for the vortex core at r m 0, the profiles f(r) for 
g = 1000, 10000 match well with the Thomas-Fermi profile 

/tfM = \J S ^2T-- Here R = ( 4 9/ 7T ) 1/4 is the Thomas-Fermi 
radius in a two dimensional trap. 

The Bogoliubov modes are solutions of the Bogoliubov- 
de Gennes equations 

~uj m v m = H G pv m + g(<j)o) 2 u m . (20) 

Here H GP = -^V 2 + \r 2 + 2 3 |0 O | 2 - M- The symmetry 
of 4>q ~ e l9 implies that 
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u(r, 9) 
v(r, 6) 



U{r) e ue+2lf > 
V(r) e tW . 



(21) 



The Bogoliubov modes can be classified into different 
multiplcts with definite angular momentum I. As only 
the modes with I = can contribute to dp{r = 0) in 
Eq.(18), from now on we focus only on the I = multi- 
plct. 

The Bogoliubov-de Gennes equations were numerically 
solved for I = using the vortex solutions </>o (Fig.l). In 
the numerical calculations the envelopes U(r) and V(r) 
were expressed as finite linear combinations of the 2D 
harmonic oscillator eigenfunctions up to certain energy 
cut-off. All the numerical quantities presented below 
were checked for convergence with the increasing energy 
cut-off. We also tested our numerical code in the limit of 
g <C 1 where some perturbative analytical results can be 
compared with numerics. The Bogoliubov modes of the 
I = multiplet can be classified as follows. 

• The anomalous mode which is a harmonic trap 
counterpart of the vortex zero mode (6) . In the trap 
the mode acquires a non-zero negative frequency 
iv a < , where u a — > in the limit of g — > oo [9] . 
It was tested in [10]. 

• The collective oscillations mode with the trap fre- 
quency u>i = 1, 

\ar r / 



n = \f + f - + rf 

. dr r 



(22) 



The frequency uj\ of the numerical counterpart of 
this exact mode was kept within 0.01 of the ex- 
act uji = 1. Contribution of the trap mode to the 
density of depleted atoms in the core can be imme- 
diately estimated as 



Mr = o)| 2 = 



(r = 0) 



1 



O(g ) . 
(23) 



This contribution does not depend on g. For large 
g it is negligible as compared to phonon modes (see 
next point). Even for a finite g this trap mode is 
a collective excitation of the center of mass of the 
atomic cloud which has nothing to do with the mo- 
tion of the vortex with respect to the condensate. 
For these two reasons we do not include the trap 
mode in dp(r = 0) even for finite g. 

Phonon modes with positive frequencies to m >i > 0. 
These modes give a nonzero contribution to dp(r = 



0) which grows with g as ^fg. The phonon contri- 
bution dominates for large g and it results in gen- 
uine fluctuation of the vortex position with respect 
to the condensate. Rather surprisingly, fluctuation 
of vortex position, if any, turns out to come from 
phonons and not from the vortex mode. 

We compared the dominating phonon contribution 
d/0 P honons(r = 0) = ^ TO >2 l w ™( r = °)| 2 to tnc quantum 
depletion in a uniform condensate without any vortex. 
If the density of the uniform condensate is equal to the 
density of our condensate near the center of the trap, 
then the quantum depletion in the uniform condensate 
is rfpuniform = g 1 ^ 2 / 4:i: 3 ^ 2 . We compare our numerical 
value of cipphonons {r = 0) to this uniform depletion in 
Fig. 2 where we plot the reduction factor 



c = 



dpphononsjr = 0) 
^Puniform 



(24) 



For large g the reduction factor e stabilizes at approxi- 
mately 0.23. We conclude that the density of depleted 
atoms inside the vortex core is a factor of 0.23 less than 
the same depletion outside the vortex. Apparently the 
centrifugal force due to vortex rotation is quite efficient 
in throwing away from the core not only the condensate 

hu^" ftlsn flip nnn-rnnrlpTiQprl atrtms 



g 

FIG. 2. e is a ratio of the density of depleted atoms inside 
the vortex core at r = to the density of depleted atoms in 
a uniform condensate without any vortex (or depletion out- 
side the vortex core). Small e means that the non-condensed 
atoms are quite efficiently expelled from the core by the cen- 
trifugal force, e stabilizes around 0.23 for large g when the 
Thomas-Fermi condensate can be very well approximated by 
a uniform condensate, at least locally close to the vortex in 
the center of the trap, e w 0.23 is the asymptotic value of the 
reduction factor in a uniform condensate. 

To summarize we find that for a large g the density of 
depleted atoms in the core is dominated by the phonon 
contribution and is approximately given by 



dp(r = 0) w 0.23 



g 1 / 2 

4^3/2 



(25) 
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It is time to compare this phonon contribution with 
the simple estimate (16) based on the truncation to the 
vortex mode. 

Mr = 0)| 2 = 



In the Thomas-Fermi limit we estimate 
2 



j_ N 

47T p 



4Li 



(r = 0) w 0.0905 1 / 2 . The ratio 



of the true to the naive density reads 

5 = ^ r = °l « 0.11. 
|«o(r = 0)| 2 



(26) 



We conclude that in the Thomas-Fermi limit the actual 
fluctuations of the vortex position are an order of magni- 
tude less than the simple estimate based on quantization 
of the effective action or equivalent truncation of the Bo- 
goliubov theory to the vortex mode. 



VI. SUMMARY 

This paper has the following logic. We start from an 
effective action for vortex dynamics in two dimensions, 
quantize it and obtain an effective wave function for the 
position of the vortex. Dispersion of this gaussian wave 
function seems to estimate quantum fluctuations in vor- 
tex position. Next we introduce the number conserving 
Bogoliubov theory and rederive the same wave function. 
This establishes a connection between the effective action 
quantization and the Bogoliubov theory: the former gives 
the same result as the latter when the latter is truncated 
to the vortex mode. Next, in the framework of the Bogoli- 
ubov theory we show that the truncation is not unique 
and that vortex position fluctuations critically depend on 
the way of truncating to the vortex mode. It turns out 
that the vortex mode in the Bogoliubov theory is not 
independent from the phonon modes so that all modes 
have to be included in order to correctly estimate quan- 
tum fluctuations. Rather surprisingly, in the full untrun- 
catcd Bogoliubov theory quantum fluctuations of vortex 
position come from phonon modes and not from the vor- 
tex mode. What is more these phonon-induccd quantum 
fluctuations are an order of magnitude less than the sim- 
ple estimate based on the truncation to the vortex mode 
or on the quantization of the effective action. We want 
to stress this result means only that the naive quantiza- 
tion of the effective action (1) docs not work properly, 
not that the effective action is incorrect. 

We also found that the density of depleted atoms inside 
a vortex core is a factor of 0.23 less than outside the core. 
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